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Optimal transport on wireless networks 
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Abstract. - We present a study of the application of a variant of a recently introduced heuristic 
algorithm for the optimization of transport routes on complex networks to the problem of finding 
the optimal routes of communication between nodes on wireless networks. Our algorithm iter- 
atively balances network traffic by minimizing the maximum node betweenness on the network. 
The variant we consider specifically accounts for the broadcast restrictions imposed by wireless 
communication by using a different betweenness measure. We compare the performance of our 
algorithm to two other known algorithms and find that our algorithm achieves the highest trans- 
port capacity both for minimum node degree geometric networks, which are directed geometric 
networks that model wireless communication networks, and for configuration model networks that 
are uncorrelated scale- free networks. 
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The study of transport on complex networks has attracted 
a great deal of interest in recent years [1-22] . One of the 
most important problems is to determine the routes be- 
tween pairs of nodes that optimize the efficiency of the 
transport. Oftentimes, the routes that are used on net- 
works are the so-called shortest-path routes, which are the 
routes with the minimum number of hops between any two 
nodes. As the volume of the transport increases, this ap- 
proach leads to congestion or jamming of highly connected 
nodes on the network called hubs. Interest has developed 
in finding the routes that allow a given network to bear the 
highest possible amount of traffic without jamming [2-7]. 
This is done in general by defining the length of a path as a 
sum of weights assigned to the links that form it, and then 
reweighting various links to spread the traffic throughout 
the network and avoid jamming at hubs. The problem of 
finding the exact optimal routing has been shown [3, 23] 
to be iVP-hard. Recently, however, we introduced a new 
heuristic routing optimization algorithm and showed [5,6] 
that it finds near-optimal routing on both random (Erdos- 
Renyi) [24] and scale-free [25] networks. Remarkably, this 
algorithm runs in only polynomial time 0(N 3 log AT). We 
also found that optimal routing preserves the small world 
character of networks [26]. 

In this paper we use a variant of our algorithm to find 
optimal routes for transport on wireless communication 



networks, which has been the subject of a number of re- 
cent papers [7-13]. Wireless networks are described by 
variants of random geometric networks. These networks 
lack the small world effect [26] that characterizes all of 
other types of networks we have previously applied our 
algorithm to. Transport on wireless networks occurs only 
along the subset of links that are bidirectional. However, 
in order to avoid broadcasting interference, every time a 
node broadcasts an information packet the nodes at the 
ends of all of its outgoing links, whether they are bidirec- 
tional or not, are prevented from broadcasting or receiving 
packets. The variant of our algorithm we consider here ac- 
counts for these broadcasting restrictions. 

In order to study the effectiveness of our optimization 
of routing on wireless networks and to understand how 
the optimization process works, we will compare results 
given by our optimal routing algorithm (OR) with those 
obtained using the shortest path routing algorithm (SP) 
and with the algorithm introduced in Ref. [7] (KR). The 
three algorithms will be applied to two different types of 
networks: minimum node degree geometric networks that, 
as stated above, are good models for wireless communica- 
tion networks [13], and configuration model networks that 
are uncorrelated scale-free networks [27]. Note that re- 
cent studies have shown that scale- free distributions of the 
node degrees are achievable [28-30] on geometric networks 
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by a community-aware manipulation of the broadcasting 
powers of the nodes. 

Minimum node degree geometric networks are con- 
structed by uniformly distributing N nodes randomly on a 
unit square and then adjusting their broadcasting powers 
until each node achieves a minimum degree k m in counting 
only bidirectional links. They differ from random geo- 
metric networks [31] by their highly peaked degree dis- 
tribution, with an average degree (k) very close to k m i n . 
Random geometric networks, by contrast, are obtained by 
connecting all pairs of nodes situated at a geometric dis- 
tance shorter than a given threshold and are characterized 
by a binomial degree distribution. To facilitate compar- 
ison with Ref. [7], the minimum degree k m i n is taken to 
be 8. Networks generated using the configuration model 
are characterized by a much broader power-law distribu- 
tion of the node degrees, p{k) oc k~^ , and by the absence 
of any correlation between the degrees of adjacent nodes. 
The node degrees are allowed to vary between a lower cut- 
off m and the square root of the number of nodes N. In 
our simulations we used m = 2 and 7 = 2.5. For both 
models we considered networks with N between 30 and 
1600. Both unconstrained routing as well as routing with 
wireless broadcasting constraints are considered for each 
case. 

Our results show that the new variant of our algorithm 
achieves a significant improvement over the one presented 
in Ref. [7] for minimum node degree geometric networks. 
However, this improvement is not as large as we achieved 
for scale-free networks. We will argue that this reduction 
in optimization efficiency is due to constraints on rerout- 
ing imposed by the non-small world nature of geometric 
networks. 

Routing on the network is assumed to be done accord- 
ing to a static protocol which prescribes the next hop(s) 
for a packet currently at node i and whose destination is 
node t. Each node is assumed to have a packet queue 
which works on a "first-in/first-out" basis. When a new 
packet is added to the network at some node or arrives at 
a new node along its path, it is appended at the end of 
the queue. Upon reaching their destination, packets are 
removed from the network. For simplicity, we assume that 
all nodes have the same processing capacity of 1 packet per 
time step (assuming they are not inhibited by broadcast- 
ing neighbors) and that new packets are inserted at every 
node at the same average rate of r packets per time step. 
This average insertion rate characterizes the load of the 
network. The destinations of the packets inserted at node 
i are chosen at random from among the other N — 1 nodes 
on the network. The algorithm can, however, be general- 
ized for nodes with different processing capacities and for 
arbitrary traffic demands. 

Given a loop-free routing table, the betweenness b\ s '^ of 
node i with respect to a source node s and a destination 
node t is defined [32] as the sum of the probabilities of 
all paths between s and t that pass through i. The total 
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Fig. 1: (Color online) Ensemble averages of the average and 
maximum betweenness as functions of network size for mini- 
mum node degree geometric networks. Lower three sets (hol- 
low black circles, red squares, and blue diamonds) represent 
{Bavg), and upper three sets (solid black circles, red squares, 
and blue diamonds) represent (B max ). 



betweenness Bi is found by summing up the contributions 
from all pairs of source and destination nodes. The prac- 
tical way [32] to compute b[ s ^ for all i and s is as follows: 
all nodes are assigned weight 1 and then the weight of ev- 
ery node along each path towards t is split evenly among 
its predecessors in the routing table on the way from t to 
s and added to the weights of the predecessors. 

The aforementioned broadcasting restrictions are equiv- 
alent to saying that every node is processing not only the 
information packets passing through itself, but also those 
passing through its incoming neighbors. This situation 
can be accounted for by using the cumulative between- 
ness 

B cum = ^B k , (1) 

fceA/i 

where the incoming neighborhood Mi is the set formed by 
node i together with all its incoming neighbors. 

The time average of the number of packets passing 
through a given node i in the course of a time step is 



rB, 
N-l' 



(2) 



while the time average of the number of packets passing 
through its incoming neighborhood is 



N - 1 



(3) 



Without broadcasting constraints, jamming of the net- 
work occurs at the critical average insertion rate r c at 
which the average number of packets processed by the bus- 
iest node reaches unity. Consequently, r c is given by [2] 
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(4) 
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Fig. 2: (Color online) Ensemble averages of the average and 
maximum cumulative betweenness as functions of network size 
for minimum node degree geometric networks. Lower three sets 
(hollow black circles, red squares, and blue diamonds) represent 
{Bavg 1 }, an d upper three sets (solid black circles, red squares, 
and blue diamonds) represent (B^™). 



Fig. 3: (Color online) Ensemble averages of the average and 
maximum betweenness as functions of network size for uncor- 
rected scale-free networks. Lower three sets (hollow black cir- 
cles, red squares, and blue diamonds) represent {B avg ), and 
upper three sets (solid black circles, red squares, and blue dia- 
monds) represent (B moI ). 



where B max is the highest betweenness of any node on the 
network. If broadcasting constraints are considered, the 
critical insertion rate is determined by the busiest incom- 
ing neighborhood and we have 

r cum = {N _ (5 ) 

Thus, to achieve optimal routing on a wireless network, 
the highest cumulative betweenness B™™ should be min- 
imized [9]. 

The application of our algorithm to ordinary networks 
has been described in Ref. [5]. For wireless networks, the 
algorithm proceeds as follows: 

1. Assign uniform weight to every bidirectional link (SP 
routing) and compute the shortest paths between all pairs 
of nodes and the betweenness of every node. 

2. Find the node io which has the highest cumulative 
betweenness B™™ and then the node with the highest 
ordinary betweenness in its incoming neighborhood Afi . 
Increase the weight of every bidirectional link that con- 
nects the latter node to other nodes by adding half the 
initial weight to it. 

3. Recompute the shortest paths and the betweennesses. 
Go back to step 2. 

To achieve the 0(N 3 log N) running time, a binary heap 
must be used in the Dijkstra algorithm [33] for the com- 
putation of the shortest paths to reduce the time required 
to sort the nodes by distance. 

The algorithm described in Ref. [7] treats bidirectional 
links as two separate directed links and uses the cumula- 
tive betweenness of the node of origin as link weight. All 
cumulative betweennesses are initially set equal to 1 for 
the purpose of routing computation and then two rounds 
of iterations are performed. In the course of each round 
all shortest path routes are computed in the order of their 



node of origin, based on the latest values of the cumula- 
tive betweennesses. These values are updated after each 
computation of the shortest paths originating from a given 
node. This algorithm runs in time 0(N 2 ). 

Throughout the paper, the network average of the be- 
tweenness Bi is denoted by B avg , while further averaging 
over an ensemble of network realizations is indicated by 
angular brackets. Fig. 1 shows the ensemble averages of 
the network average and maximum betweenness, (B avg ) 
and (B max ) respectively, as functions of the number of 
nodes N for minimum node degree geometric networks. 
Results are presented for shortest path routing (SP), for 
the optimal routing given by our algorithm (OR), and 
for the routing algorithm described in Ref. [7] (KR). The 
results for the average cumulative betweennesses {B™™) 
and (-B™J™) for minimum node degree geometric networks 
are shown in fig. 2. For small networks of up to approxi- 
mately 100 nodes the scaling of (B™™) and (B™™) with 
network size is somewhat anomalous, but for larger net- 
works both quantities scale with network size according 
to power laws. The exponents of the power laws were 
computed by fitting data corresponding to N between 200 
and 1600. The values of the exponents are given in ta- 
ble , with the quoted errors being 2a estimates. It is 
apparent from figs. 1 and 2 and from table that our al- 
gorithm achieves the lowest maximum betweenness for a 
given network size. Both optimization algorithms obtain a 
significant improvement over shortest path routing. This 
is unlike the case of random or scale- free networks, where 
our algorithm brings a significant improvement in trans- 
port capacity when compared to any other optimization 
algorithm. 

There are two topological reasons for this less significant 
difference. The first reason is the lack of a small world 
effect [26]. In geometric networks with node-to-node com- 
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SP 


OR 


KR 


{Bavg) 

{Bmax) 
1 jDcum\ 

/ r>cum\ 
\ 7nax/ 


1.488 ±0.003 
1.874 ±0.023 
1.459 ±0.004 
1.759 ±0.005 


1.519 ±0.006 
1.585 ±0.019 
1.492 ±0.002 
1.495 ±0.005 


1.509 ±0.003 
1.619 ±0.013 
1.477 ±0.002 
1.568 ±0.007 



Table 1: Exponents of the {B avg ), {B max ), (-B£™), and 
{Bmax) power-law scaling with network size N for minimum 
node degree geometric networks. 



munication ranges much smaller than the physical size of 
the network any shortest paths follow approximately the 
geometric shortest path and pass on average through a 
number of nodes proportional to the square root of the 
geometric distance between source and destination. On 
virtually all other types of networks of practical interest 
the average number of hops along the path increases with 
network size slower than logarithmically. This absence of 
shortcuts causes network traffic to be quite evenly spread 
even in the case of SP routing and reduces the likelihood 
of finding alternative paths that lower the maximum be- 
tweenncss. Computation of (L avg ) shows that, even in the 
case of the minimum node degree model where the broad- 
casting powers of the nodes arc not equal, the average 
number of hops along the path scales with network size 
approximately as y/~N. The lack of a small world effect 
could be remedied by adding a few long range connections 
between randomly chosen pairs of nodes [34, 35] . This 
can be done in practice by scattering a few special nodes 
equipped with long range unidirectional antennas. The 
second reason for the less significant difference between 
the two algorithms is the highly peaked distribution of 
the node degrees which is a characteristic of the minimum 
node degree model and also contributes to the uniformity 
of traffic spreading. 



Next, we computed (B a 



(B„ 



[BZT), and 



(Bmax) f° r uncorrected scale- free networks generated us- 
ing the configuration model using all three types of rout- 
ing. Scale-free networks are characterized by an average 
number of hops which increases with the number of nodes 
slower than logarithmically. Results are shown in fig. 3 for 
(B avg ) and (B max ) and in fig. 4 for (B£J™) and (B™£). 
The exponents of the power laws for (B max ) and (B^™) 
are given in table . The quantities (B avg ) and (B%£™) 
vary in this case approximately as YlogY [6] and fit- 
ting them with a power law doesn't make sense from a 
theoretical point of view. The factor of improvement pro- 
vided by our routing optimization algorithm is even more 
significant than in the case of the minimum node degree 
geometric network model and is comparable for both wire- 
less and ordinary networks. This may become particularly 
useful if the geometric network topology of wireless net- 
works is somehow altered (for example by introducing a 
community-aware distribution of the broadcasting powers) 
to resemble the topology of a small world network. 



{Bmax} 
/ T?cum\ 
\ D raaxl 



SP 

1.626 ±0.011 
1.799 ±0.010 



OR 
1.184 ±0.012 
1.480 ±0.013 



KR 
1.308 ±0.010 
1.617 ±0.008 



Table 2: Exponents of the {B max ) and (B%££) power-law scal- 
ing with network size N for uncorrelated scale- free networks. 



In summary, we have introduced a new algorithm for 
transport optimization on wireless networks and compared 
its performance with another recently introduced routing 
optimization algorithm. The effectiveness of transport op- 
timization on wireless networks was compared to results 
obtained for ordinary networks with no broadcasting con- 
straints. We found that our algorithm performs better in 
all cases studied, more significantly so in the case of scale- 
free networks. The less significant difference in the case 
of minimum node degree geometric networks is due to the 
geometric network topology (which results in a lack of a 
small-world effect) and to the narrowly peaked distribu- 
tion of node degrees. 
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